Mon. Not. R. Astron. Soc. 000, 1-13 (2006) Printed 5 Febraary 2008 (MN I^TbX style file v2.2) 



Radio and infrared recombination studies of the southern massive 
star-forming region G333.6— 0.2 

Takuya Fujiyoshi/'^* Craig H. Smith,^* James L. Caswell/* Toby J. T. Moore,^* 
Stuart L. Lumsden,^* David K. Aitken''* and Patrick F. Roche^* 

^Subaru Telescope, National Astronomical Obsen'atory of Japan, National Institutes of Natural Sciences, 650 North A'ohoku Place, Hilo, HI 96720, USA 
^Department of Physics, University College, University of New South Wales, Australian Defence Force Academy, Canberra ACT 2600, Australia 
■^Electro Optic Systems Pty Limited, 55A Monaro Street, Queanbeyan NSW 2620, Australia 
^Australia Telescope National Facility, CSIRO, P. O. Box 76, Epping NSW 2 12 1, Australia 

^Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Egerton Wharf, Birkenhead CH41 ILD 
^Department of Physics and Astronomy, University of Leeds, Leeds LS2 9JT 
^Department of Physical Sciences, University of Hertfordshire, Hatfield, Herts ALIO 9AB 
^Astrophysics, Department of Physics, Oxford University, Keble Road, Oxford 0X1 3RH 



J> ■ To be published in MNRAS 



ABSTRACT 

We present high spatial resolution radio and near-infrared hydrogen recombination line obser- 
vations of the southern massive star-forming region G333.6— 0.2. The 3.4-cm continuum peak 
is found slightly offset from the infrared source. The H90q! spectra show for the first time a 
double peak profile at some positions. The complex velocity structure may be accounted for 
by champagne outflows, which may also explain the offset between the radio and infrared 
sources. The 2.17-|j.mBr7 image and H90q! map are combined to construct an extinction map 
which shows a trend probably set by the blister nature of the H 11 region. The total number 
of Lyman continuum photons in the central 50-arcsec is estimated to be equivalent to that 
emitted by up to 19 07V stars. 
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1 INTRODUCTION 

G333.6— 0.2 is a southern compact Hll region which was ini- 
tially detected as one of the brightest radio sources in a 5-GHz 
Galactic plane survey (Goss & Shaver 1970). Radio recombina- 
tion lines (RRLs) H90a and He90a were observed by McGee, 
Newton & Batchelor (1975), and McGee & Newton (1981) ex- 
amined H76a and He76a RRLs. Both papers reported a Local 
Thermal Equilibrium (LTE) electron temperature (T*) of slightly 
less than 7000 K and a ratio of singly ionised helium to singly 
ionised hydrogen [Y^ = ji^'^+y] of less than the solar value (i.e., 
< 10 per cent). The H II region is almost totally obscured in the 
visible {Ay = 18.5, Landini et al. 1984) but is bright at infrared 
(IR) wavelengths, and has been studied intensively in this region 
of the spectrum (e.g. Aitken, Griffiths & Jones 1977; Hyland et al. 
1980; Fujiyoshi et al. 1998, 2001, 2005, hereafter Papers 1, 2, and 
3, respectively). In Paper 1, it was found that the Hll region must 
be excited by a cluster of O and B stars, rather than by a single star. 
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Indeed, near-IR (NIR) broadband imaging revealed multiple point 
sources in G333.6— 0.2 (Paper 3). Point sources with significant 
H — K colour excess were found preferentially in extended emis- 
sion regions and it was suggested that these red objects were phys- 
ically associated with the H II region. The powerful nature of this 
nebula was unveiled by mid-IR (MIR) imaging polarimetry (Pa- 
per 2) which showed bent magnetic field lines possibly shaped by 
the expansion of the ionised gas. 

At radio wavelengths the most intense emission from H II re- 
gions is caused by interactions between unbound charged particles, 
better known as free-free radiation or bremsstrahhing (braking ra- 
diation). Since the particles are free and their energy states are not 
quantised the radiation resulting from changes in their kinetic en- 
ergy is continuous. Free-free emission is dominated by interactions 
of free electrons with the most abundant element, singly ionised 
hydrogen. Observations of the radio free-free emission (where ex- 
tinction is usually negligible) can be used to predict the hydrogen 
recombination line fluxes at IR wavelengths, which in turn can be 
compared to the observed values to estimate extinction at the line 
wavelengths. 

In this paper, we present radio synthesis observations at 
3.4 cm, and the 2.17-fi.m Br7 hydrogen recombination line imag- 
ing of G333.6— 0.2. High spatial resolution data provided by radio 
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interferometry can be compared directly with those at NIR wave- 
lengths to infer spatial variation of extinction across the face of the 
H II region. The H90a RRL, observed both at high spatial and spec- 
tral resolution, is used as a probe of the physical conditions in the 
emitting gas. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Radio observations 

The radio observations were made on 1996 October 8 at the Aus- 
tralia Telescope Compact Array (ATCA) in Narrabri, Australia. The 
ATCA is a synthesis telescope which consists of six 22-m antennae 
on a 6-km east-west track. One set of measurements was taken with 
the baseline configuration 6A. The longest baseline of 6 km was 
chosen to give the final map the highest possible spatial resolution 
available with the ATCA. 

G333.6— 0.2 was observed at a frequency centred between the 
H90a (rest frequency 8.872 GHz) and He90a (8.876 GHz) RRLs. 
Each of two orthogonal polarisations was recorded, using a corre- 
lator configuration of 512 channels covering a 16-MHz bandwidth. 
The primary calibrator (1934-638; 2.8 Jy at 8.8 GHz with an un- 
certainty of up to a few per cent) was observed for 20 minutes dur- 
ing the run for the bandpass correction and the flux calibration. A 
secondary calibrator (1619—680) was observed every 20 minutes 
for 5 minutes throughout the observation period for position and 
phase calibration. 

Data reduction was carried out using the software package 
MIRIAD. The continuum and H90a RRL maps were constructed 
using natural weighting in the conventional Fast Fourier Trans- 
formation. The maps were then deconvolved using the CLEAN 
routine in MIRIAD. Spectral cubes were also made following the 
same method; however, in order to improve the signal-to-noise ra- 
tio (S/N), 4 consecutive channels were binned together resulting in 
an effective velocity resolution of ~ 4 km s^^. The cleaned maps 
were finally restored to intensity maps with a 1.6-arcsec diameter 
circular beam, which therefore is the effective spatial resolution of 
the maps presented here. 



2.1.1 Missing flux from extended structures 

McGee et al. (1975) observed G333.6-0.2 as part of their H90a 
RRL survey of bright southern radio sources. Using a 2.5-arcmin 
beam they measured the 3.4-cm continuum flux density of 80 ± 
17 Jy. However, in our continuum map the integrated total flux den- 
sity in an area ~ 2 x 2 arcmin^ was found to be only ~ 13 Jy. 

Lack of short baselines in interferometry leads to the loss of 
sensitivity to extended structures. The ATCA 6A configuration has 
a shortest baseline of 337 m, and thus at 3.4 cm is not respon- 
sive to structure smoothly distributed on a scale of ~ 30 arcsec 
or larger. The target is a compact H II region whose IR counterpart 
has a core-halo geometry (Aitken et al. 1977) in which most flux is 
concentrated in the central region. We have assumed that the miss- 
ing extended emission underlying the compact H II region takes the 
shape of a 2-D Gaussian. Its peak value and full widths at half max- 
imum (FWHMs) are adjusted so that it would contain the missing 
flux and satisfy the level of extinction measured by Landini et al. 
(1984) (see Section 3.2.1 for details). 



Table 1. Summary of NIR observations of G333.6-0.2. 



Date 


Filter 


Ao 


AA 


Integration 






iyxm) 




(sec/position) 


11 Aug 1995 


Br7 


2.17 


1 % 


2 X 10.0 




H2(2-1)S(1) 


2.25 


1 % 


1 X 10.0 



2.2 NIR observations 

The NIR observations of G333.6— 0.2 were carried out on 1995 
August 11. All measurements were obtained using the IR imag- 
ing spectrometer IRIS at the //36 Cassegrain focus of the 3.9-m 
Anglo-Australian Telescope (AAT), Siding Spring, Australia. The 
detector element of IRIS is a Rockwell 128 x 128 HgCdTe NIC- 
M0S2 array, which is housed in a closed-cycle compressed he- 
lium refrigeration unit. In the imaging mode, wavelength selection 
is provided by discrete filters such as 'standard' NIR broadband 
(J, H, K). A more complete description of the instrument can be 
found in Allen et al. (1993). 

Table 1 summarises the NIR observations of G333.6— 0.2. IN- 
TERM mode was used to give a pixel size of ~ 0.25 x 0.25 arcsec^. 
The object was dithered in the instrument's field of view, which 
produced the final image size of ~ 50 x 50 arcsec'^. The 2.25- 
\im H2(2— 1) S(I) image was used for subtracting the contin- 
uum from the 2. 17-|i-m Br7 image. Moneti & Moorwood (1989) 
found that G333.6— 0.2 showed only a weak H2(l— 0) S(l) emis- 
sion in the central region [/bit = 3.5 X 10"^** W cm and 
/(i_o) = 1-3 X 10^^" W cm^^, both measured in a 6 x 6-arcsec^ 
beam]. In conditions characteristic of H II regions, the (1—0) mode 
is always much stronger than the (2—1) (e.g. Burton 1992). We 
therefore assumed that the (2—1) image contained only the con- 
tinuum emission. The calibration was achieved by observing the 
standard star SA94-251 {K = 8.305, Carter & Meadows 1995). 
Seeing on the night was quite good and observations of standard 
stars showed FWHMs of about 0.6 arcsec. 



3 RESULTS AND DISCUSSIONS 
3.1 Radio data 

3.1.1 The 3. 4-cm continuum 

Figure 1 shows the 3.4-cm continuum map. The overall similarity 
to the IR images (cf. Figure 7. See also Papers 1 & 3) is quite re- 
markable. It is dominated by a north-south elongated main peak in 
the centre with steeper gradient to the east. Hyland et al. (1980) 
suggested that this steep gradient may be caused by extinction. In- 
deed, the extinction map (see Section 3.2.1) shows a slight increase 
of extinction in this direction. An arm-like extension, which was 
described as the secondary peak in Papers 1 & 3, is seen about 
5 arcsec to the east of the main peak. The gap between the main 
and secondary peaks is also present. As the gap is still present at 
3.4 cm, it is likely that the main and secondary peaks are separate 
ionised regions and the gap represents a narrow non-ionized region 
in between, rather than an apparent structure simply caused by a 
heavy foreground extinction. It was suggested in Paper 3 that both 
main and secondary peaks house exciting sources. 

The most significant difference between the IR and radio con- 
tinuum images is that the radio continuum main peak is found off- 
set slightly (~ 2 arcsec) from the MIR main peak position. The 
coordinates of the radio main peak are RA = 16''22™09^5 and 
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Figure 1. The 3.4-cm continuum grey-scale map of G333.6— 0.2. Contours 
are drawn at -0.050 (broken line), 0.00050, 0.10, 0.15, 0.20, 0.30, 0.40, 
0.50, 0.60, 0.70, and 0.79 Jy beam"! . Axis offsets are from the MIR main 
peak position [RA = 16'"22™09?6 and Dec = -50°06'00" (J2000)]. 
The circular 1 .6-arcsec diameter synthesised beam is also shown at the bot- 
tom left hand corner. 



Dec = -50°05'58" (J2000). The size of the source, determined 
by fitting a 2-D Gaussian profile, is 4.9 arcsec in RA and 6.9 arcsec 
in Dec, or about 0.086 pc ~ 1.8 x 10'' au at 3.0 kpc (an assumed 
distance to G333.6— 0.2, Colgan et al. 1993). There is a ridge run- 
ning to the south (and slightly east) from the main continuum peak 
and in fact the MIR main peak position coincides with the peak on 
this ridge. 



3.1.2 TheH90aRRL 

The H90q RRL map is presented in Figure 2. It shows the spa- 
tial variation of the integral over velocity of the recombination line. 
The distribution of the RRL emission is similar to that of the 3.4-cm 
continuum. The size of the H90a peak emission, which is coinci- 
dent with the peak of the 3.4-cm continuum emission, measured 
again by fitting a 2-D Gaussian profile, is 4.1 arcsec in RA and 
5.7 arcsec in Dec, which is only slightly smaller than the size of the 
continuum peak. There is also a ridge running to the south of the 
main peak; however, the peak on the ridge is slightly more south- 
east (by about 1 arcsec south and 0.5 arcsec east) than the ridge 
peak position in the 3.4-cm continuum map. 

Figure 3 shows the total H90a RRL spectrum in the central 
~ 2 X 2 arcmin^ . Attempts were made to fit either a single or a dou- 
ble Gaussian profile to the spectrum and it was found that a combi- 
nation of two Gaussian functions gave a better fit (the sum of fitting 
residuals was smaller). As it will be shown later, G333.6— 0.2 ex- 
hibits a complex velocity structure with some positions displaying 
a double peak spectrum. Such a structure has never been observed 
with large-beam, single-dish radio studies of the H II region and this 
emphasises the importance of high spatial resolution investigations. 

The parameters derived from the fit are summarised in Table 2. 
All velocities refer to the Local Standard of Rest (LSR) frame. 



RA offset (arcsec) 

Figure 2. The H90a grey-scale map of G333.6— 0.2. Contours are drawn 
at -0.25 (broken line), 0.025, 0.25, 0.65, 1.3, 1.7, 2.1, 2.5, 3.0, 3.5, 4.0, 4.3 
Jy-km s~^ beam~^. The coordinates of the origin are as in Figure 1. 



-0.2 




Figure 3. The total H90q spectrum in the central ~ 2 X 2 arcmin . It 
shows the data points (dots), a multiple (double) Gaussian fit to the data 
(solid line), and the residual (i.e., the difference between the actual data and 
the fit; dashed line). 

The weak broad feature centred at around —270 km s~^ is the 
HI 13/3 RRL. Additional weak emission may be present at about 
— 150 km s"'^ and could be the He90Q RRL, although it would 
normally be expected at a bluer velocity [~ —175 km s~^, cf. Fig- 
ure l(o), McGee et al. (1975)]. It should be noted, however, that 
both lines presumably also exhibit double peak spectra and there- 
fore are likely to be blended. No attempt was made to fit parameters 
to these RRLs. 

Table 2 also lists the single-dish H90q measurements of 
McGee et al. (1975). The overall velocity of the ionized gas in 
G333.6— 0.2 was found by three independent single-dish observa- 
tions at different frequencies to be ~ —46 km ^T^ . Apart from 
McGee et al. (1975), who found its velocity to be —47.1 km 
McGee & Newton (1981) determined -45.3 km s"^ from the 14.7- 
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Table 2. The Gaussian fit parameters to the total H90q: spectrum in the 
central ^2x2 arcmin^ (see Figure 3). 



of the true electron temperature, Te, using (Brown 1987) 

1 0.87 



Si 






Reference 


(mjy) 


(kms-i) 


(km s-i) 




433 ±117 


-73.4 ±3.1 


26.3 ±4.5 


This study 


798 ± 53 


-45.6 ± 2.8 


36.1 ±4.2 


This study 




-47.1 ± 0.2 


43.9 ±0.3 


McGee et al. (1975) 



GHz H76a RRL observation with a 2.3-arcmin beam, and Caswell 
& Haynes (1987) obtained —46 km s^^ by measuring the 5-GHz 
H109a and HIIOq? RRLs with a ~ 4-arcmin beam. The redder 
(~ —46 km s~^) of the two velocity components we found above 
is therefore likely to be associated with emission from extended 
structures, most of which escaped detection in the current study 
(see Section 2.1.1). The bluer component at ~ —70 km s~^, de- 
tected for the first time with high spatial resolution observations 
presented here, must then be originating from compact structures. 

Figure 4 shows average spectra in 1.75 x 1.75-arcsec^ boxes 
in the central region superimposed on the 3.4-cm continuum map. 
The boxes are positioned so that they would be centred on the main 
peak and the peak on the ridge. Positions 34, 35, and 36 are centred 
on each of the three radio peaks in the south-west extension (see 
Figures 1 & 2). Figure 5 shows the H90a RRL spectra fitted with 
Gaussian profiles. Again, we determined whether each spectrum 
contains a single or a double line by fitting either a single or a dou- 
ble Gaussian profile and selecting the one that resulted in a smaller 
sum of residuals. Parameters found by the fit are summarised in 
Table 3. 

Table 4 lists central velocities summarised in Table 3 in a con- 
figuration closer to that shown in Figure 4 so as to illustrate the 
spatial trend of the velocity fields. It clearly shows that the single 
and double peak regions are segregated. At and to the north-west 
of the peak of the H90a RRL emission (Position 13) the single line 
has the mean central velocity of ^ —53 km s^^. Double peaks al- 
most surround the above single peak region with the mean central 
velocities at ~ —65 km s^^ and ~ —34 km s^^. 



3.1.3 Temperatures 

The LTE electron temperature is given by (Roelfsema & Goss 
1992) 



T* 



6943 ■ p 



Si ■ AVg 1 + F^ 



where f is the frequency in GHz, Sc and 5"; are continuum and line 
flux densities in Jy beam^^, respectively, and AVg is the Gaussian 
line width in km s^^. Shaver et al. (1983) found to be nearly 
constant in the Galaxy, and we adopted this value (0.074). 

The brightness temperature, Tb, is given by (e.g. Wood & 
Churchwell 1989) 



5clO- 



where c is the speed of light, k is the Boltzmann's constant, and S7 
is the beam size. Ti, is related to Te by 

Ti, = Te(l-e-"=), 

where Tc is the continuum optical depth. 

may be corrected for non-LTE effects to derive an estimate 



1 - 



where bn and I3„ are the LTE departure coefficients, which depend 
on Te and Tie (electron density), and can be computed using the 
code devised by Brocklehurst & Salem (1977). At each position, 
we calculated , Tt, and Tc. We then created a grid of 6„ and j3„ 
for a range of Te (2000 K ^ Te ^ 12000 K, ATe = 100 K) and 
rie (1 X 10^ cm"^ ^ Tie ^ 1 X 10^ cm-^ Aue = 1000 cm"^). 
Applying these correction factors, we selected a combination of 6„ 
and f3„ (hence n^) that gave the closest corrected Te to the input 
Te. These values are listed in Table 3 and in Table 4 again in a 
configuration similar to that shown in Figure 4. 

The mean values of Te and Ue found in this way are Te = 
5700 ± 1300 K and Ue = 407000 ± 226000 cm'^*. The result 
implies the fractional mean-square Te variation [t^ = (AT/T)^, 
Ferland (2001)] for G333.6-0.2 is ~ 0.052. While is thought to 
be ~ 0, since a gas in photoionisation equilibrium is nearly isother- 
mal (Kingdon & Ferland 1995), relatively high values of have 
been found in some nebulae (e.g. 0.032 in M8, Esteban et al. 1999; 
0.028 in Orion, O'dell, Peimbert & Peimbert 2003). However, it 
should be noted that electron temperatures derived from measure- 
ments of RRLs, when corrected for the spatial variation of Y^, 
tend to show smaller t^. For example, Roelfsema, Goss & Mallik 
( 1992) detected a large variation of Y^ in W3A, and compensating 
for this fluctuation found the mean LTE electron temperature to be 
T* = 7500 ± 750 K, or t'^ = 0.01. We did not derive Y+ but 
assumed a constant value (0.074, Shaver et al. 1983). 



3.2 The Br7 image 

Figure 6 shows the 2.17-^m Br7 plus the underlying continuum. 
There are at least 1 1 point sources recognisable by eye in the field 
imaged, and positions of 6 of these sources coincide with those of 
NIR sources found common in all the broadband JHK' frames 
(Paper 3). It is interesting to note that, as a simple analysis of rela- 
tive positions amongst the point sources reveals, it is not the bright- 
est point source in Figure 6 that coincides with the MIR main peak 
(Paper 1), but the one immediately below [i.e., the source at offset 
(0,0)]. 

Figure 7 shows the 2.17-nm Br7 image after the continuum 
subtraction, which has been Gaussian smoothed to a slightly lower 
effective spatial resolution of ~ 0.8 arcsec. The morphology seen 
in Figure 7 is strikingly similar to that in the MIR continuum im- 
age (see Figure 2, Paper 1). G333.6— 0.2 has been found to have a 
blister geometry viewed face-on (Hyland et al. 1980), and the sim- 
ilar morphology exhibited by the ionised gas and warm dust may 
be accounted for if the MIR emitting dust grains are located behind 
the H II region on the front surface of the photo-dissociation region 
(Paper 1). It is also possible that the gas and dust are well mixed 
but it is generally thought that dust grains are sublimated at about 
1500 K. 

Table 5 lists the measured Br7 fluxes as a function of the aper- 
ture sizes and compares the results of the present study to the previ- 
ously published measurements. Each measurement was made with 
the aperture centred on the main Br7 peak [offset (0,0)] and the 
overall uncertainty is estimated to be about 5 per cent. The results 
obtained are consistent with those from the previous studies, espe- 
cially with the measurements of Wynn- Williams et al. (1978), and 
imply that the continuum subtraction process was reasonably suc- 
cessful. 



Recombination studies of G333.6—0.2 



Table 3. The Gaussian fit parameters to the H90a RRL from the grid of 1.75 X 1.75-arcsec^ boxes (see Figure 5). The position numbers refer 
to those indicated in Figure 4. Si is in mJy beam^^ and rie in X 10^ cm^^. 



No. Si err Vlsr err AVg err T* T5 Tc Te err rie bn Pn 

(km/s) (km/s) (K) (K) (K) 



1 


9.3 


4.1 


-65.7 


7.4 


35.8 


9.8 


4861 


856 


0.19 


5000 


4000 


3.00 


0.997977 


-0.372641 




5.4 


2.0 


-29.4 


18.8 


46.1 


30.1 


















2 


14.0 


1.2 


-68.1 


0.9 


20.3 


2.2 


3166 


488 


0.17 


3200 


800 


1.65 


0.997272 


-0.184034 




6.1 


0.7 


-31.4 


3.7 


42.3 


9.1 


















3 


16.9 


1.2 


-70.8 


1.0 


24.1 


2.4 


4477 


883 


0.22 


4500 


1200 


4.58 


0.998264 


-0.069765 




6.3 


1.0 


-24.9 


3.4 


40.1 


9.6 


















4 


8.0 


0.9 


-59.4 


2.5 


21.1 


4.2 


3093 


209 


0.07 


3100 


1600 


1.65 


0.997268 


-0.149526 




3.6 


1.0 


-38.7 


5.0 


19.5 


9.3 


















5 


23.7 


1.3 


-60.2 


1.9 


29.9 


4.1 


3834 


935 


0.27 


3900 


2100 


2.51 


0.997812 


-0.160766 




6.6 


2.7 


-34.2 


4.5 


19.5 


8.4 


















6 


30.7 


9.4 


-69.7 


2.7 


27.7 


4.2 


4350 


2269 


0.59 


5100 


2100 


1.90 


0.997535 


-0.691416 




21.8 


3.8 


-41.3 


8.2 


41.5 


11.5 


















7 


26.6 


1.9 


-68.7 


1.3 


27.1 


3.0 


4349 


1773 


0.47 


4700 


900 


2.44 


0.997793 


-0.404504 




18.1 


1.7 


-26.6 


2.2 


36.0 


5.6 


















8 


7.2 


1.5 


-71.7 


3.3 


28.9 


8.6 


4980 


1101 


0.25 


5000 


1700 


7.71 


0.998473 


-0.049835 




17.1 


1.5 


-23.0 


1.3 


30.5 


3.4 


















9 


3.2 


1.1 


-77.5 


7.1 


35.9 


20.5 


4780 


718 


0.16 


4800 


2900 


5.66 


0.998362 


-0.078471 




12.6 


1.2 


-22.8 


1.5 


30.3 


4.3 


















10 


8.6 


0.3 


-56.8 


0.5 


24.0 


1.4 


6063 


422 


0.07 


6200 


2300 


2.51 


0.997843 


-0.800377 




2.0 


0.5 


-37.5 


1.1 


8.7 


2.7 


















11 


19.0 


0.9 


-56.5 


0.7 


33.6 


1.8 


4279 


810 


0.21 


4300 


300 


3.95 


0.998179 


-0.070785 


12 


45.9 


1.4 


-55.6 


0.7 


43.6 


1.6 


4761 


2869 


0.70 


5700 


200 


2.51 


0.997835 


-0.664589 


13 


71.3 


2.2 


-57.6 


0.7 


49.9 


1.8 


4169 


4381 


1.09 


6600 


300 


1.62 


0.997404 


-1.284620 


14 


34.6 


1.8 


-71.1 


1.0 


27.0 


2.2 


6011 


3716 


0.68 


7500 


1000 


4.09 


0.998172 


-0.853680 




27.7 


1.4 


-30.9 


1.5 


37.9 


4.0 


















15 


17.5 


1.9 


-72.0 


1.5 


23.9 


3.7 


7540 


3461 


0.47 


9200 


1700 


6.18 


0.998316 


-1.098160 




29.5 


1.7 


-29.4 


1.1 


34.0 


3.0 


















16 


5.7 


1.2 


-78.5 


3.5 


30.4 


9.7 


5914 


1654 


0.29 


6600 


2200 


2.37 


0.997801 


-0.949274 




22.7 


1.1 
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Figure 4. A grid of average spectra in 1.75 X 1.75-arcsec^ boxes superimposed on the 3.4-cm continuum map of G333.6— 0.2. Contours are drawn from 
0.1 Jy beam^^ in an 0.1 Jy beam^^ interval. 



Table 5. Measured Br7 fluxes (in X 10 ^® W cm as a function of beam 
size (arcsec). 



Beam diameter 


Br7 flux 


Br7 flux 




[this study] 


[previous studies] 


6.0 


3.2 


2.8" 


6.5 


3.6 


3.5'' 


30 


14 


7.3" 


32 


15 


12'' 



" Landini et al. (1984); Wynn-WilHams et al. (1978) 
3.2.1 Extinction map 

The theoretical Br7 flux (in W cm~^) is related to the radio flux 
density Sv (Jy) measured at frequency v (GHz) by 



8.90 X asv^-^V^^ SyTl 



where as is the total recombination coefficient to the n = 2 level 



(the 'on-the-spot' approximation), and F = 



is the number 



ratio of Lyman continuum to Br7 photons. V can be computed for 
a range of Te and using tables of Storey & Hummer (1995). 

The observed flux density is related to the theoretical value in 
the usual way, i.e., 

j-,obs _ ritheo -T2.17 

where superscripts obs and theo represent observed and theoretical 
fluxes, respectively, and T2.17 is the optical depth at the wavelength 



of the Br7 emission. The extinction at a given wavelength A is re- 
lated to the optical depth by 

A\ = 2.5rA logio e. 

Finally, assuming the extinction at 2.17 nm, Ai.\-j ~ Ak, the ex- 
tinction found at 2.17 (im can be extrapolated to the visual extinc- 
tion at 0.55 H-m using the NIR extinction law derived by Martin & 
Whittet(1990) 



Ai 



: O.OMv 



For Te, the nearest to the mean value derived earlier (5700 K) 
in the table of Storey & Hummer (1995), 5000 K was used. The 
mean value of obtained at the same time showed quite a large 
standard deviation (He = 407000 ± 226000 cm"^); however, the 
2. 17-)i.m Br7 emission line is not very sensitive to electron densi- 
ties in the range usually found in H II regions. Populations of lev- 
els are more susceptible to density in high n states, such as those 
found at radio frequencies, as incoming free electrons have a higher 
probability of perturbing the bound electrons into a different orbit. 
Whereas low n transitions [such as Br7 (n = 7—4)] are affected 
less by collisions. We nonetheless used a mean value in the density 
range ne = 10'' — 10^ cm~'^ in the calculation of the theoretical 
Br7 flux. 

Some positions exhibit a moderate optical thickness at 3.4 cm 
(see Table 3). 5*1, at each position has therefore been multiplied by 
acorrection factor, [tc/(1 — e^^'=)], which varied from 1.03 atone 
end of the region mapped with a reasonable 5'/A'^ (Position 17) to 
1.64 at the peak (Position 13). However, the correction can only be 
made in regions where Tc is available, i.e., only in the 36 1.75 x 
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Figure 5. The H90a RRL spectra averaged in 1.75 X 1.75-arcsec^ boxes (see Figure 4). Each window shows data points (dots), a Gaussian profile fit to the 
data (solid line), and the number at the top left hand corner refers to the position number in Figure 4. 



1.75-arcsec boxes. This will confine the extent of the extinction 
map to the central ~ 10 x 10-arcsec . 

Landini et al. (1984) measured Av — 18.5 mag in a 6- 
arcsec aperture centred on the 7v -band peak. As discussed in Sec- 
tion 2.1.1, radiation from extended structures escaped detection in 
our current radio interferometric observations. We assumed that 
the distribution of the missing flux underlying the central compact 
structure is a 2-D Gaussian. We then adjusted its peak value and 
FWHMs so that it would contain the amount of the missing flux 
and contribute sufficient emission in the centre to yield the level 
of the Landini et al.'s extinction measurement. This was achieved 
by adding a 2-D Gaussian with a peak value of 0.4 Jy beam"^ 
and FWHMs of 20.5 arcsec; however, we note that recently re- 
leased Spitzer GLIMPSE (Benjamin et al. 2003) images (IRAC 3.6, 
4.5, 5.8, & 8.0 (xm) of the G333.6— 0.2 region show a wide spread 
(~ 2 X 2 arcmin^), fairly uniform faint emission around more in- 
tense (in fact, saturated) central object. Therefore, it is quite pos- 
sible that there is a large scale pedestal-like extended component 
to the missing flux, in addition to the 2-D Gaussian structure. If 



so, the peak value and FWHMs we used for the Gaussian missing 
flux component are likely to be their respective upper and lower 
limits. Consequently, our map is more suited to demonstrate the 
spatial trend in the extinction rather than to show the exact extinc- 
tion values. It is also noted that, although altering the peak value 
and FWHMs of the 2-D Gaussian missing flux component changed 
the finer details of the resultant extinction map, the overall spatial 
trend remained virtually the same in the central region. Even adding 
a uniform missing flux produced a very similar trend. 

Figure 7 shows the extinction map (contours) superimposed 
on the Br7 image of G333.6— 0.2. Lower extinction values are 
found near the Br7 peak [offset (0,0)] and in the western exten- 
sion [around offset (—6,-1)], i.e., near and around the ionising 
sources. Point sources with the brightest apparent K magnitudes 
and the reddest H — K colours were found at these positions, and 
the large H — K colour excess of those objects most probably in- 
dicates the physical association of the point sources with the H II 
region (Paper 3). This trend seems to fit with the blister nature of 
G333.6-0.2. 
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Table 4. Spatial distribution of central velocities of the H90o RRLs found by the Gaussian profile fit to a spectium averaged in a 1 . 75 X 1 . 75-arcsec'^ box 
(see Figure 5). Each box lists, from the top to bottom, the position number as in Figure 4, the central velocity (velocities if a double peak profile is found at 
the position) in km s^^, and non-LTE corrected electron temperature (in K) and electron density (in cm^'*). 







4 

-59.4 ±2.5 
-38.7 ±5.0 
3100 
165000 


10 

-56.8 ±0.5 
-37.5 ±1.1 

6200 
251000 


17 

-56.3±11.2 
-42.8±13.6 

6600 
485000 










2 

-68.1 ±0.9 
-31.4±3.7 

3200 
165000 


5 

-60.2 ±1.9 
-34.2 ±4.5 
3900 
251000 


w 

-56.5 ±0.7 

4300 
395000 


18 

-53.4 ±0.6 

5000 
714000 


24 

-52.9 ±0.8 

5200 
371000 






1 

-65.7 ±7.4 
-29.4±18.8 
5000 
300000 


3 

-70.8 ±1.0 
-24.9 ±3.4 
4500 
458000 


(j 

-69.7 ±2.7 
-41.3 ±8.2 
5100 
190000 


12 

-55.6 ±0.7 

5700 
251000 


19 

-54.0 ±0.7 

6400 
476000 


25 

-51.3 ±0.7 

6400 
899000 


29 

-54.1 ±2.0 

5900 
96000 








7 

-68.7 ±1.3 
-26.6 ±2.2 
4700 
244000 


13 

-57.6 ±0.7 

6600 
162000 


20 

-56.7 ±0.7 

6400 
319000 


26 

-52.3 ±0.8 

7500 
254000 


30 

-48.5 ±2.0 

5200 
364000 
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-71.7 ±3.3 
-23.0 ± 1.3 

5000 
771000 


14 

-71.1 ± 1.0 
-30.9 ±1.5 
7500 
409000 


21 

-66.8 ±3.8 
-37.3 ±4.1 
7000 
191000 


27 

-57.2 ±2.1 
-34.6 ±1.7 
5800 
862000 


31 

-47.0 ± 1.4 

5900 
48000 


33 

-58.4±33.8 
-38.4±11.7 
4100 
252000 
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-77.5 ±7.1 
-22.8 ±1.5 

4800 
566000 


15 

-72.0 ±1.5 
-29.4±1.1 

9200 
618000 


22 

-71.3 ±2.1 
-36.7 ±1.2 

8800 
880000 


28 

-62.6 ±3.0 
-35.6 ±1.0 
6700 
630000 


32 

-61.7±5.0 
-36.4 ±1.4 
5700 
331000 
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-78.5 ±3.5 
-30.0 ±0.8 
6600 
237000 


23 

-55.5 ±6.0 
-38.0 ±0.7 
5500 
607000 


34 

-61.2 ±4.5 
-35.3 ±1.4 
6000 
308000 




36 

-64.7 ±5.0 
-41.7 ±3.5 
6000 
556000 


Position 
Vlsr 
Vlsr 












35 

-67.5 ±3.9 
-41.3 ±3.4 
4700 
580000 





A considerable number of H II regions have been found to 
have a blister geometry in which the H 11 regions are located at the 
edges of molecular clouds (Habing & Israel 1979). The blister ge- 
ometry is thought to be the result of the champagne phase in which 
H II regions expand preferentially toward directions of decreasing 
density (e.g. Yorke 1986). The blister H II regions are therefore ion- 
isation bounded on the molecular cloud side and density bounded 
on the side of outward champagne flow (Yorke, Tenorio-Tagle & 
Bodenheimer 1983). The density of ionized gas decreases approx- 
imately exponentially away from the cloud and a large fraction of 
ionizing photons can escape from the density bounded side. That 
is, that side can be relatively optically thin. Hyland et al. (1980) 
found that G333.6— 0.2 is a blister H II region viewed face-on. The 
density bounded side then is the front face and the low extinction 
values found around the ionizing sources are perhaps natural con- 
sequences of this geometry. 

Taken at face value, the presence of red objects in lower ex- 
tinction regions may seem contradictory; however, it was shown 
in Paper 3 that, in the colour-colour diagram space, for example, 
these red sources are found significantly away from the average in- 



terstellar extinction curves. It was also demonstrated that emission 
from hot (> 600 K) dust grains, which are likely to be located in 
circumstellar discs, is considerably more important in the A'-band 
than at H, resulting in the large H ~ K colours of these objects. 
In other words, the red colours are due not to the extinction but to 
the intrinsic infrared excess, and so it is not implausible for these 
sources to be found in lower extinction regions. 

There is a slight increase in extinction between the main and 
secondary peaks. In a clumpy medium, an expansion front (weak 
ionisation front with its associated shock) of an H II region passes 
around the more dense clumps of neutral material and compresses 
them (Yorke 1986), which results in embedded neutral globules and 
dark lanes sometimes seen in optical images of emission nebulae. 
As mentioned earlier, there is likely to be an ionizing source both in 
the main and secondary peaks and the gap between the two peaks 
observed at NIR and MIR wavelengths could be a pressure com- 
pressed dense region that consequently shows a high extinction, 
rather than a foreground extinction lane. The fact that this gap is 
also observed at radio frequencies (see Figures I & 2) may support 
this suggestion. 
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high extinction. The position coincides with a radio pealc (see Sec- 
tion 3.3) and the high extinction may be associated with this radio 
source. 
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Figure 6. A grey-scale image of the 2.17-^m Br7 plus the continuum. Co- 
ordinates of the origin are as in Figure 1 . 
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Figure 7. The 2.17-iam Br7 image of G333.6— 0.2. Grey scale: Same as 
Figure 6 but after the continuum subtraction and a Gaussian smoothing (see 
Section 3.2). Contours: Extinction map (see Section 3.2.1). Contours are 
drawn at Ay = 4 mag interval from 12 mag (near the origin) to 36 mag 
[near (-3,4)]. 



The highest extinction is located near offset (—3,4). It was 
suggested in Paper 2 that there is likely to be an extra overlapping 
of dense, cold material containing silicates at 5 arcsec north of the 
main MIR peak, since the spectrum taken at the position showed a 
deep silicate absorption. The beam size used for the 5'.' ON position 
was 5.6 arcsec in diameter and most likely included this region of 



3.2.2 The Lyman continuum photons 

The theoretical Br7 flux can also be expressed as 



Br7 



Br7 



he 



47rD2 ABr7 



where D is the distance to the object (in cm), h is the Planck's 
constant, and ABr7 is the wavelength of the Br7 emission. And 
recalling that F = 



Nbt 



NLy 



he 



Again by knowing S,^, Ue and Te, and using an appropriate table 
from Storey & Hummer (1995), the number of Lyman continuum 
photons can be estimated. 

In constructing the extinction map, we have assumed a 2-D 
Gaussian-shaped missing flux distribution. This also means that the 
missing flux is concentrated and that ~ 100 per cent of the total 
flux at 3.4 cm (80 Jy, McGee et al. 1975) comes from the inner 
region. However, such a high concentration of flux in the centre is 
most probably the upper limit as there appears to be a wider spread 
of faint emission further out (see above). The number of Lyman 
continuum photons in the inner ^ 50-arcsec is thus estimated to be 
< 9.5 X 10*^ s-'K 

Rubin et al. (1994) suggested that the effective temperature 
(Teff) of star(s) exciting G333.6-0.2 is ~ 36000 K which, accord- 
ing to the stellar spectral-type scale of Martins, Schaerer & Hillier 
(2005), is approximately characteristic of an 07V star. The number 
of the Lyman continuum photons emitted by an 07V star is esti- 
mated to be ~ 5 X 10** s~^ (Martins et al. 2005). Therefore, if 
all the exciting stars in G333.6— 0.2 were of the same spectral type, 
there would be up to 19 07V stars. In Paper 1, it was found that 
about a dozen 08V stars {T^s = 36000 K) is required to radiate 
the ionising photon flux necessary to produce the observed [Ne II] 
central dip.^ 



3.3 Velocity structure (revisited) 

The double peak of an RRL observed in an H II region could be 
interpreted as: 

(a) expansion of shell-like structure, 

(b) rotation, 

(c) two H II regions at different velocities, or, 

(d) outflows/jets. 

We discuss each of these in turn, followed by a short section (e), 
describing details of a toy model for our preferred option, (d). 

(a) expansion of shell-like structure. There is no obvious shell- 
like morphology seen in the regions where double peaks are ob- 
served. The mean velocity difference between the two lines is ~ 
30 km s^^. For a shell-like structure to be unresolved (< 1.6 arcsec 



^ Note that the stellar parameter calibration of Martins et al. (2005) takes 
into account effects of line-blanketing which shifts (lowers) the effective 
temperature for each stellar spectral type. However, for a given stellar ef- 
fective temperature, other stellar parameters (e.g. luminosity, etc.) are es- 
sentially unchanged. 
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or < 0.023 pc) and to have an expansion velocity of ^ 30 km s^^ 
would require a dynamical age of only a few hundred years. This 
seems rather unlikely. Also, the stellar-wind blown cavity scenario 
for the [Ne II] central dip was ruled out in Paper 1, since there were 
no signs of such a cavity in the continuum image at a nearby wave- 
length. 

(b) rotation. This appears unlikely since there is no systematic 
shift of central velocities across the region to suggest a rotation. 

(c) two H II regions at different velocities. Combined with the 
multiple point sources seen in the NIR images (see Figure 6 and 
Paper 3), it could be that the double peak profiles of the H90q! RRL 
in G333.6— 0.2 represent the existence of several compact Hll re- 
gions at different velocities. In the JHK' images (Paper 3), all of 
which were taken at a similar spatial resolution to that of the ra- 
dio maps presented here, at least four point sources were identified 
in the central region of G333.6— 0.2, two of which are found in 
the main peak (cf. NM and SM in Figure 7, Paper 3), one in the 
secondary peak to the east of the main, and one in the south-west 
extension. Double peak regions do indeed coincide with the MIR 
main peak and the south-west extension. At the highest spatial reso- 
lution so far attained 0.6 arcsec), the 2. 17-|i.m image (Figure 6) 
reveals even more star-like point sources in the central region and 
there appears to be three sources located very close together near 
the position of the MIR main peak. 

Figure 8 shows intensity distribution in each velocity plane ('^ 
4 km s~^ wide) in the range ~ —100 to km s^^. This sequence 
of velocity maps also reveals the multiple nature of the exciting 
source in the central region of G333.6— 0.2. The main radio peak 
dominates throughout the velocity range because of its large line 
width 50 km s^^ FWHM). The ridge to the south of the main 
peak is visible up to the —62 km s~^ frame but is absent till another 
peak appears further south at about —50 km s~^. The south-west 
extension seems to contain at least three peaks (see also Figures 1 & 
2). As mentioned earlier, there is at least one point source identified 
in the NIR images in the south-west extension (see Figure 6 and 
Paper 3). 

The shape of the radio main peak is rather complex. There ap- 
pears to be another peak close (about 2 arcsec north-west) to the 
main peak (see the —58 km s^^ frame). The extinction map (Fig- 
ure 7) shows a large increase in the extinction value around this 
position, which may imply that obscuring dust grains are associ- 
ated with this peak. The main peak and this nearby local maximum 
do not have obvious IR counterparts and could be younger than the 
other exciting sources in the H II region, which in turn may also 
explain the presence of heavy obscuring dust grains. If this suppo- 
sition of the difference in age were correct then this alone would 
lead to the conclusion of multiple exciting sources, and hence the 
presence of several compact H II regions, in G333.6— 0.2. 

If flat Galactic rotation curves, with a constant circular ve- 
locity Q = 220 km s^^ and the Sun's distance from the Galac- 
tic centre i?o ~ 8.5 kpc, were adopted, the difference in velocity 
of ~ 30 km s~^ between two peaks in the double peak regions 
would imply a separation of the H II regions along the line of sight 
(los) of ~ 1.6 kpc. While possible, it seems rather unlikely that 
several unassociated compact H II regions just happen to coincide 
along the los to produce the velocity structure observed. Of course, 
G333.6— 0.2 is an active star-forming complex and adopting the 
Galactic rotation curves to interpret the velocity structure may not 
be appropriate. 

Another interesting feature that is difficult to interpret with 
the los coincidence scenario for the complex velocity structure 
presents itself when the single and double peak spectra are com- 
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Figure 9. Gaussian profiles at Positions 13, 14, and 15 (see also Figure 5). 

pared closely. Figure 9 shows the best-fit Gaussian profiles at Po- 
sitions 13, 14, and 15. It is remarkable that the overall shape of 
the spectra matches almost exactly and the double peak appears 
more or less symmetric with the single peak between the two peaks. 
Again it seems rather unsatisfactory to explain such features as a 
chance coincidence of compact H II regions along the los. 

The spatial resolution of the radio maps (1.6 arcsec), despite 
being the highest so far obtained at these frequencies and for the 
first time revealing the double peak profile of the H90a RRL, is 
more than twice the highest resolution of the IR images (~ 0.6 arc- 
sec). To examine more closely the relationship between the IR point 
sources and the velocity structure, radio interferometric observa- 
tions at even higher spatial resolution, at least that comparable to 
the IR images, are required. IR imaging both at high spectral and 
spatial resolution is perhaps an alternative. 

(d) outflow or jets. There are no obvious structures that suggest 
outflows or jets in the double peak regions. However, it is possible 
that outflows/jets are inclined in the los so that such structures are 
not immediately obvious, or are unresolved by the 1.6-arcsec beam. 
Yorke et al. (1983) produced theoretical radio continuum maps of 
H II regions in the champagne phase, which occurs when the ion- 
isation front of an H II region encounters a region of strong den- 
sity gradients, such as the boundary of a molecular cloud. They 
found that the champagne outflow (or preferential expansion of an 
H II region along negative density gradients) can attain high veloc- 
ity (> 10 km s^^) and in the case of a slab-like geometry it can 
be bipolar. They also found that in general the maximum of radio 
emission, which is proportional to (e.g. Rubin 1968), does not 
coincide with the ionizing source. This may explain the offset be- 
tween the radio and IR sources. 

The champagne outflow scenario may also explain the broad 
line width of the single peak. The outflows caused by the prefer- 
ential expansion of the H II region into less dense regions may not 
be well-coUimated. In fact, the velocity structure depicted in the 
models of Yorke et al. (1983) shows a fan-like radial outpouring 
of ionized gas from the mouth of the outflow. If an RRL were ob- 
served from such a structure its spectrum would probably show a 
large line width due to a range of los velocities present in the radial 
velocity field. 

It was found in Paper 2 that the intrinsic magnetic field lines 
follow well the flux distribution north-west of the MIR main peak. 
This morphology is most likely to be due to the presence of a clump 
of dense material containing silicates (indicated by the shape of the 
spectrum taken at 5 arcsec north of the MIR main peak). There is 
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Figure 8. The H90a RRL maps in each velocity plane (~ 4 km s ^ wide). The velocity plane is indicated at the top left hand corner of each plot. The contours 
are drawn at a 5 per cent interval from to 100 per cent of the maximum flux density in each plane. The maximum values in each plane are 1.4, 1.8, 2.2, 2.3, 
3.8,5.0,5.8, 8.1,8.0,8.3,8.4, 9.2,7.5,6.9, 6.0,5.3,4.0,3.7,3.2,2.8, 2.5, 1.5, 1.9, and 0.82 xlO'^ Jy beam"! from -100 to -2.9 km s"!. 
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Figure 10. A 'cartoon' representation (a side view) of a possible situation 
in G333.6— 0.2. denotes the IR exciting source(s). Two different lines 
of sight (1 and 2) are also indicated, together with a probable RRL profile 
along each los. Arrows are drawn only to illustrate likely direction of the 
champagne flows. See text for details. 



4 CONCLUSIONS 

The radio interferometric observations and NIR imaging of the 
southern massive star-forming region G333.6— 0.2 are presented. 
The morphologies of the high spatial resolution (1.6-arcsec) 3.4- 
cm continuum and H90a RRL maps are remarkably similar to each 
other and to the IR images. However, the radio main peak neither 
coincides with the IR main peak nor has any obvious IR counter- 
parts. The H90q RRL spectra show a complex structure with a dou- 
ble peak profile at some positions observed for the first time. The 
complex velocity structure may be explained by champagne out- 
flows, which may also explain the offset between the radio and IR 
sources. The 3.4-cm radio continuum map is combined with the 
2.17-fi.m Br7 image in constructing the visual extinction map. The 
blister geometry of the H II region probably causes the trend seen 
in the extinction map. The number of Lyman continuum photons in 
the central 50-arcsec is estimated to be equivalent to that emitted 
by < 19 07V stars. 



a local visual extinction maximum about 4 arcsec north and 3 arc- 
sec west of the IR main peak (Figure 7), which may be associated 
with the radio peak ~ 2 arcsec north and ~ 3 arcsec west of the 
radio main peak seen in the —58 km s~^ plane (see Figure 8). The 
expanding Hll region may have carried the frozen-in, threading 
magnetic field with its expansion. If such structure then encoun- 
tered the aforementioned clump of dust and compressed the mag- 
netic field against it the intrinsic polarization pattern may result. 
Near the boundary, where the expansion hit the clump, the radio 
continuum peak may have been produced by the high density. 

(e) a toy model. In view of the ability of the outflow model 
to account for the general appearance, we present (Figure 10) a 
more elaborate toy model of G333.6— 0.2, which can explain all 
of the key features. The blister geometry, that may account for the 
trend found in the extinction map (Figure 7), may not be exactly 
face-on but somewhat tilted south-east. Towards the north of the 
exciting (IR) source(s), as we look through a range of los velocities 
in the fanned-out champagne outflow in the relatively high den- 
sity regions near the molecular cloud wall, we observe the strong, 
broad single RRL (los I, Figure 10). We obtain double RRL spec- 
tra towards the south-east, as faster flows are found along the wide 
opening of the blister 'bowl', whilst flows in the intermediate ve- 
locity range are generally fainter (as the region is more diffuse) and 
relatively perpendicular with respect to the los (los 2, Figure 10). 
Champagne flow velocities are expected to be greater further away 
from the molecular cloud (Bedijn & Tenorio-Tagle 1981) and we 
are possibly observing this phenomenon as there seems to be an 
apparent trend that the blue peak shifts to the bluer velocity (and 
the red redder) the more south we look (see Table 4). 

Another interesting feature, which may be explained by this 
picture, is that in the north the blue peak is stronger, whereas the red 
peak becomes brighter in the south (see Figures 4, 5, and 9). This 
could be due to the presence of the dense clump. The density near 
the front side of the molecular cloud wall is enhanced by the dense 
clump and therefore, near the clump, we find strongest emission 
from the flows that are directed towards us; south of the clump, 
emission from the back wall (away from us) dominates. Note that, 
whilst the strength of the blue peak changes dramatically depending 
on the exact position, the red peak brightness remains more or less 
constant. 
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